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Abstract 



The problem of extrapolation and interpolation of asymptotic series is consid- 
ered. Several new variants of improving the accuracy of the self-similar approxi- 
mants are suggested. The methods are illustrated by examples typical of chemical 
physics, when one is interested in finding the equation of state for a strongly in- 
teracting system. A special attention is payed to the study of the basic properties 
of fluctuating fluid membranes. It is shown that these properties can be well de- 
scribed by means of the method of self-similar approximants. For this purpose, the 
method has been generalized in order to give accurate predictions at infinity for a 
function, whose behavior is known only at the region of its variable close to zero. 
The obtained results for fluctuating fluid membranes are in good agreement with 
the known numerical data. 
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1 Introduction 



Asymptotic series are ubiquitous, arising in practically all realistic problems that do not 
allow for exact solutions but require the use of some kind of perturbation theory. The 
latter assumes that there exists a small parameter, such that the observable quantities 
of interest can be represented as series in powers of this asymptotically small parameter. 
However, in applications, this parameter is not negligibly small, but takes finite values. 
In the majority of cases, the perturbative series are divergent for the parameter finite 
values, corresponding to realistic systems. This is why the problem of extrapolating 
asymptotic series is of great importance. The problem poses the question what are the 
most accurate ways of extending the validity of the series, obtained for asymptotically 
small parameters, to the finite values of the latter. Moreover, in some cases it is necessary 
to extend perturbative results to the extreme limit, where the parameter tends to infinity. 

The most popular techniques of extrapolating asymptotic series are based on the Pade 
approximants [1]. However, these have several well known shortcomings. First of all, they 
are not uniquely defined. For a series of a given order, there is a whole table of many 
Pade approximants, but there is no general recipe that would advise which of them should 
be preferred. The conclusion could be made if all these approximants from the table 
would be close to each other. Then their dispersion would define their accuracy. This, 
unfortunately, is a very rare case, since a rather standard situation is when the results of 
the approximants from the table are widely scattered, making it difficult to judge on their 
accuracy. Even worse, among the approximants from the table, there very often appear 
those that contain unreasonable poles. Then the accuracy of such approximants, strictly 
speaking, is not defined at all. One can, of course, neglect the approximants with the 
poles, calling them outliers. But, as is clear, this is a quite subjective procedure, since it 
may happen that the arising poles could have meaning, as is the usual case for critical 
phenomena. Arbitrarily rejecting some of the approximants, while keeping others, renders 
the whole procedure not well defined. Then one is not able to ascribe any accuracy to 
the obtained results. The most difficult situation is when one needs to extrapolate the 
perturbative series to the infinite value of the expansion parameter. In such a case, the 
Pade approximants cannot be used at all. In order to be specific, showing that this 
necessity of extrapolating the parameter to infinity does happen in realistic cases, we can 
mention the problem of calculating the pressure of fiuctuating fiuid membranes. 

Different types of membranes are rather frequent structures in chemical and biological 
systems [2-6]. The membrane thermal fluctuations between two hard walls are often 
described by field theory. An important class of membranes arc fluid membranes, whose 
constituent molecules are able to move within them. The fluctuations are controlled by 
their bending rigidity. 

When applying fleld theory to the description of membranes, one encounters the fol- 
lowing problem. First, to proceed in calculations, one replaces the hard walls by a smooth 
potential of a finite stiffness, which, in dimcnsionless units, can be denoted as g. Then, 
one is able to proceed by invoking perturbation theory in powers of g. However, to return 
back to the sought case of rigid walls, one needs to set 51 ^ oo. Thus, the problem arises 
how from an expansion in powers of small g ^ one could extract information on the 
quantities of interest for large g — )■ 00? The standard resummation techniques, such as 
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Borcl or Fade [6] could help in extending asymptotic series in small — >■ to finite values 
of g, but these techniques are not applicable for the limit of (7 — )■ cxd. 

The problem of extending a function /((?), which is known only for asymptotically 
small — )■ 0, to the whole region of g, including the limit g ^ 00, can be solved by the 
optimized perturbation theory, advanced in Ref. [7]. This theory has been successfully 
applied to a variety of problems, as can be inferred from the review-type articles [8,9] and 
references therein. The pivotal idea of the optimized perturbation theory [7] is to introduce 
in the calculational process control functions defined by optimization conditions. As a 
result, a function f{g), known only for small g ^ 0, can be extrapolated to the whole 
region of g, including the limit g ^ 00. There exist three main ways of introducing 
control functions. One way is to include in the initial approximation trial parameters 
that are transformed, by means of the optimization conditions, into control functions at 
each step of perturbation theory. Another way is to introduce control functions in the 
process of accomplishing a perturbative or iterative scheme, for instance, by defining the 
cutoffs of integrals or introducing regularization masses, which are then to be transformed 
into control functions. The third way is to derive, first, an asymptotic series in powers 
of a small parameter or variable and then to reorganize the derived series by means of 
a change of variables, with control functions included in this variable transformation. 
Various examples of introducing control functions can be found in literature [10-30] (see 
also the review articles [8,9]). 

The optimized perturbation theory has been applied to fluid membranes in several 
papers. The most accurate results have been obtained by Kastening [31], whose calcula- 
tions are based on the sixth-order perturbation theory, with introducing control functions 
by means of the Kleinert change of variables [32]. This method requires rather heavy 
numerical calculations. Also, the Kastening result [31] shghtly deviates from the Monte 
Carlo [33] simulations for this problem. 

The aim of the present paper is to reconsider the general problem of extrapolating 
asymptotic series in order to find an accurate and simple way of accomplishing such an 
extrapolation. We also consider the case when the behavior at infinity is known, which 
then becomes the problem of an accurate interpolation. We pay a special attention to 
developing simple methods for obtaining the limit /(oo) at ^ — > 00 from the asymptotic 
expansion oi f{g) ai g ^ 0. To illustrate the methods, we apply them to several problems, 
which yield the series, whose mathematical structure is typical for many real situations 
in chemical physics. At the end, we apply these methods to considering the problem of 
fluid fluctuating membranes. 

The method we aim at developing is based on the self-similar approximation theory 
[34-41] in the variant involving the self-similar root approximants [42-44] and self-similar 
factor approximants [45-49]. This approach has been applied to variety of problems, 
providing high accuracy and at the same time being quite simple [42-49]. It was shown 
to be essentially more accurate than the use of Fade approximants [42-52]. However, 
in some cases, when the self-similar approximants were directly applied for solving the 
problems involving the limit of g — )■ 00, the results were not satisfactory. Our aim now 
is to generalize the method of self-similar approximants so that it would provide good 
accuracy for the limiting value of /(oo) at g ^ 00. Then we apply the method to 
solving several problems requiring the construction of the equations of state for strongly 
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interacting systems. Among them, we study the problem of fluctuating fluid membranes 
and demonstrate that the newly developed methods provide good accuracy by comparing 
our results with the known numerical calculations. 



2 Factor and root approximants 

First, we need to recall the method of self-similar factor and root approximants [42-52], 
which we are going to improve. Suppose, we are interested in the behavior of a real 
function f{g) of a real variable g. And let us assume that this function is defined by 
so complicated equations that the sole thing we are able to find is the property of the 
function at asymptotically small g, where 

m ^ fk{g) {g ^ 0) (1) 

is approximated by the series 

A(5) = /o(^)Ea.5", (2) 

n=0 

in which /o(5') is a form that cannot be expanded in powers of g. Without the loss of 
generality, we can set Oq = 1, since any oq not equal to one can be incorporated into fo{g). 

Series (2), which are valid for g ^ 0, can be extrapolated to the region > by means 
of the self-similar factor approximants [45-49] having the form 

m9)-Mg)Ui^+A^9r , (3) 

i=l 

where 

^ U/2, A: = 2,4,... 

\ {k + l)/2, k = 3,5,... 

The parameters Ai and rii are defined by the re- expansion procedure, when the k-th order 
approximant (3) is expanded in powers of g up to the k-th order as 

rk{9)^fo{9)j:a:g\ (5) 

n=0 

where a* = a* ({Aj}, {rij}). Then expansions (2) and (5) are compared, with equating the 
same-order terms 

<({A,},K})=a„. (6) 

This way is also called the accuracy-through-ordcr procedure. 

In order to give an explicit representation of Eqs. (6), it is convenient to equate the 
logarithms 

lnf;{g)c^]nMg) {g ^ 0) , (7) 
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which, taking into account form (3), yields 



Nk k 

^n,ln(l + A^)~ln^a^^-. (8) 

i=l m=0 



Expanding here 

oo / i\m— 1 



ln(l + A^)=E (^i9r^ 



m=l 

we come to the equations 



m 



J2n,A^ = B^ {n=l,2,...,k) , (9) 



1=1 

in which 



(-l)n-i r/" / " \ 

^n = )-^lim — In E^-^™ • (10) 



(n - 1)! fl->o dg"" 



\m=0 



The system of equations (9) contains k equations. When k is even, system (9) defines 
all k/2 parameters Ai and k/2 parameters n^. When k is odd, then there are k + 1 
parameters, {k + l)/2 parameters Ai, and {k + l)/2 parameters rii. Then the system (9) 
is complemented by the condition Ai = 1 following from the scaling arguments [47,48]. 
Thus, the re-expansion procedure completely defines all parameters of the self-similar 
approximant (3). It may happen, though it is a rather rare case, that Eqs. (9), for some 
order k, do not have sohitions. Then one just needs to proceed to the higher orders of the 
series. But it is important to stress that for each given order k the factor approximants 
are uniquely defined. So, when Eqs. (9) possess solutions, these solutions are unique. 

The described above method allows for the extrapolation of a series for a small g ^ 
to the finite values oi g. In some cases, there can exist additional information on the 
behavior of the function at asymptotically large g, 

fig) ^ fp{9) ^ oo) , (11) 

so that 

p 

fpig) = E bjg'"' (ttj > , (12) 

i=i 

with the powers aj in the descending order. In that case, we have the problem of inter- 
polation between small g ^ and large g ^ oo. Suppose that fo{g) at large g behaves 
as 

fo{g) ~ Ag'^ (g^oo). (13) 

Then, in order that the factor approximant (3) would satisfy the limiting form (12), we 
need to set 

AY[Ai' ^bi a + ^rii^ ai . (14) 

1=1 i=l 

This type of the interpolating factor approximant, employing k terms from the small- 
variable expansion and one limiting term from the large- variable behavior, will be denoted 
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as fk+i{g)- In the following sections, we develop alternative methods of interpolation, 
improving the accuracy of the factor approximants 

The problem of interpolation can also be conveniently solved by involving the self- 
similar root approximants [42-44], having the form 

Kig) = Mg) (. . . ((1 + A,gr + A^g'Y + . . . + A.g^Y" . (15) 

The parameters Ai and rij are defined by the large- variable expansion (12). Again, it is 
important to emphasize that this definition is unique [9]. If, instead, we try to define 
the parameters of a /c-order root approximant by the accuracy-through-order procedure, 
expanding Eq. (15) in powers of g and equating the resulting expansion with Eq. (2), 
then we confront the problem of nonuniqueness of solutions for the sought parameters 
[50]. In the following sections, we shall suggest a way of solving this problem. It is worth 
noting that the regions of small g and large g can be easily interchanged by the change 
of the variable g tol/g. 

3 Problem of self-similar interpolation 

One of the well known difficulties in dealing with asymptotic series occurs when the 
number of their terms is small, which does not allow one to construct higher order ap- 
proximants. In the present section, we suggest a way of overcoming this difficulty. The 
method, we advance, reminds the learning algorithms used in statistical learning [53]. 
The idea is as follows. Suppose we have k terms ak of the small- variable expansion and 
the limiting form of the large-variable behavior. Interpolating from the right to left, that 
is, considering the variable l/^f, we construct the corresponding self-similar approximant, 
say, the root approximant Rl{g)- Then we expand the latter in powers of g up to the 
{k + l)-order, obtaining an additional term al_^i. Using the new expansion, we define the 
approximant Rl_^^{g). Expanding this up to the (/c + 2)-order, we find the {k + 2)-order 
term 0^+2- Then, we construct the approximant Rk+2{g)i ^^d so on. Thus, each approx- 
imant defines the higher-order term of the small-gf expansion. Of course, this procedure 
can work only when the sought function pertains to the class of monotonic functions and 
the interpolation problem is considered. Below, we illustrate the method by examples 
where one is interested in finding the equations of state. 

3.1 Frohlich optical polaron 

Let us consider the problem of the optical polaron, being interested in finding its energy 
e{g) as a function of the coupling parameter g. It is common to employ the dimensionless 
notations for these quantities, which we use in what follows. The sm.dX\-g expansion and 
the large-^f limit can be found in the review article [54]. For the small-^f expansion, one 
has 

e{g) ~ aig + 02/ + 03/ (^ ^ 0) , (16) 

with the coefficients 

Qi = -1 , 02 = -1.591962 X 10-2 , 03 = -0.806070 x 10"^ , 
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While the large- behavior is given by the Miyake limit 



e{g)c^Bg^ + 0{l) ^ oo) , (17) 

where B = -0.108513. 

Following the procedure described above, we derive the coefficients 

al = -5.014168 x 10"^ , 4 = -3.312472 x 10"^ . 

Using these and interpolating from the right to left, with the variable l/g"^, we construct 
the root approximant 



V I ( I ( 1 + -^y +-^Y+ -^y +^]"' , (18) 



where 



Al = 64.163254, A2 = 7.001856 x 10^ 
As = 7.026125 x 10^ A4 = 5.201706 x 10^ 



3 5 7 1 

and the root approximant 



with the coefficients 

Al = 68.38553, A2 = 7.742967 x 10^ ^3 = 8.213401 x 10^ 

^4 = 7.313112 x 10^ A5 = 4.417553 x 10^ 

3 5 7 9 1 

The accuracy of these approximants can be checked by comparing them with the results 
of the Monte Carlo simulations [54] accomphshed for the region oi g E [1, 15]. In all this 
region, the approximant (19) has the percentage error less than 1%. The maximal error 
of —1% occurs at g = 10, where this error is comparable with the Feynman variational 
calculations [55]. But for all other values of g in the considered interval, the accuracy 
of approximant (19) is better than the Feynman results. Comparing the accuracy of the 
approximants R^g), Rl{g), and -^5(5'), we observe numerical convergence. For instance, 
the maximal percentage error of -^3(5') is 1.5% at 51 = 10. 
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3.2 One-dimensional Bose gas 



Let us now consider the ground-state energy e{g) of the Lieb-Liniger model [56] as a 
function of the couphng parameter g, again using dimensionless units. The weak-couphng 
expansion can be written [57,58] as 

e{g) ~ ^ + a^g''/^ + a^g^ + a^g''^^ , (20) 

with the coefficients 



03 = -0.424413, 04 = 0.065352, 05 = -0.017201. 
For strong couphng, we have the Tonks- Girardeau hmit 

e{g) ^'^ + o(^^ ^ 00) . (21) 

Following the procedure, described at the beginning of this section, we find 

og = 5.153629 X 10"^ 

Then we construct the root approximants of different orders, interpolating from the right 
to left, with the variable 1/g, and compare their accuracy with numerical data [58]. The 
approximant Rl{g) has the maximal, with respect to the whole range of g E [0, 00), error 
of 3.4% at g = 6. The maximal error of Rl{g) is 1.75% at g = 10. The best approximant 
is obtained using the coefficients 03 and 04 of expansion (20) and the three terms of the 
strong coupling limit that can be written [57,58] as 

e{g)c^^(l - -+^^] (g^oo). (22) 



9 9^ 



The corresponding root approximant is 



where 

Ai = 8.126984, A2 = 37.345427, A3 = 164.914098, 

A4 = 388.171278, A5 = 385.382911, 

3 5 7 9 1 

rii - , "-2 = ^ ' 6 ' ^ 8 ' ^ 5 ■ 

This aproximant (23) provides a very high accuracy in the whole range of (? G [0, 00), 
having the maximal error of only 0.023% at g = 6. Therefore, expression (23) can be 
employed as an analytical representation for the ground-state energy of the Lieb-Liniger 
gas. 
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3.3 Diluted Fermi gas 



Let us now turn to defining the ground-state energy of the spin 1/2 Fermi gas with attrac- 
tive interactions, corresponding to the negative scattering length Ug. The ground-state 
energy e(g) can be written [59] as an asymptotic expansion in powers of the dimensionless 
parameter g = Ikpcisl, where A;^ is the Fermi wave vector, 

e{g) ~ oo + aig + Osg'^ + aag^ + a^g'^ . (24) 

Here — > and 

3 1 

oo = TTT' "^1 = ^2 = 0.055661, 

10 Stt 

ag = -0.00914, 04 = -0.018604. 
In the unitary hmit [60,61], when g oo, numerical calculations [62] yield 

Jim e(g) = 0.132 . (25) 

We construct the self-similar approximants of different orders and compare their ac- 
curacy with Monte Carlo simulations [63]. The factor approximant f^+iig) turns out to be 
analogous to the diagonal [2/2] Fade approximant [64]. However, the root approximant 
R^{g) is essentially more accurate. The factor approximant f4^i{g) displays the same 
accuracy as R^{g). The maximal percentage error of the latter two approximants in the 
interval of g E [0, 5] is only about 0.2%. 



4 Problem of self-similar extrapolation 

The problem of extrapolation of asymptotic series is much more difficult than that of 
their interpolation. In the latter case, the large- variable limit is given, while in the former 
case, this limit is not known. And often, it is exactly the limiting behavior at large 
variable, which is of the most interest. In the present section, we suggest new variants 
of constructing the self-similar approximants in the extrapolation problem and illustrate 
these methods by the model whose mathematical structure is typical of the variety of 
physical and chemical systems. 



4.1 Iterated root approximants 

The extrapolation of asymptotic series can be done by means of the self-similar factor 
approximants. But, as is mentioned in Sec. 2, if we try to accomplish the extrapolation 
by using the root approximants, we encounter the problem of nonuniqueness of defining 
their parameters by the accuracy-through-order procedure. To overcome this problem, 
we suggest to use the iteration method, by keeping the fixed lower-order parameters when 
constructing the higher-order approximants. Then all parameters of the root approxi- 
mants can be uniquely defined. 

To be concrete, let us consider the anharmonic-oscillator model with the Hamiltonian 
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where x e (—00, 00) and the dimensionless couphng parameter is positive, being in the 
region g e [0, 00). The weak-couphng expansion for the ground-state energy [65] reads as 

e{g) ~ ao + aig + 025^^ + 0351^ + 045^^ + a^g^ + a^g^ + oyfif^ , (27) 

with the coefficients 

1 3 

oo ^1 ^ 4' = -2.625, 

aa = 20.8125, 04 = -241.2890625, ag = 3580.98046875, 

as = -63982.8134766, 07 = 1329733.72705. 

Our aim is to extrapolate the weak-couphng expansion (27), vahd for asymptoticaly small 
— > 0, to the region of finite values of g. And we shall pay a special attention to 
the behavior of the extrapolated energy at large — >■ 00, comparing it with the known 
asymptotic form 

e{g) ~ 0.6679865^/''^ (^ ^ 00) . (28) 

The iteration method for constructing the uniquely defined root approximants is eluci- 
dated in the following forms extrapolating the asymptotic expansion (27) to finite values 
of g. We start with the lowest approximant 

Rl{9)-\{l + Mgr , (29) 

in which 

Ai = 8.5, ni = 0.176 . 

The next-order approximant is 

Rl{9) = \{{.^ + Mgr + A^g'Y' , (30) 

with the same and 

^2 = 227.719, n2 = — = 2.771, = 0.064. 

The next approximant 

^6(^) = ^ (((1 + A^)- + A25T' + ' (31) 

contains the same Ai, n2, and A21 with 

^3 = 3.827 X 10^ 714 = — = 2.001, ng = 0.032. 

77,5 

We shall compare the strong-coupling behavior of these root approximants with the exact 
limiting form (28) and with the factor approximant 

rM = ^(1 + B,gr^{l + B^gr^l + B,gr^ , (32) 
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in which 

Bi = 26.74018, B2 = 12.46882, B3 = 3.83804, 
mi = 1.80165 X 10~^ ms = 0.05473, mg = 0.20047. 
The strong-couphng behavior of the root approximants (29) to (31) is 

Rlig) ^ , Rlig) - ^ {A,g^T ' ^^l^) - ^ (^2^^)"" • (33) 

While the factor approximant (32), as gi — > 00, gives 

f;{g) ~ ^ Bp g"'^+^^+^^ . (34) 

Substituting here the corresponding values of the parameters yields 

Rl{g) ~ 0.728698^°-^^^ , Rl{g) ~ 0. 70769 lc/°-^''^ , 

Rl{g) ~ 0.707814y°-i92 , f;{g) ~ 0.756157^°-25^ . (35) 

Comparing these expressions with the exact asymptotic form (28), we see that the am- 
plitudes of the root approximants provide slightly better extrapolation than the factor 
approximant, however this difference is not essential, all amplitudes being defined with 
an error of about 10%. In many cases, the most important quantity that is required 
to be found from the extrapolation procedure is the power of g in the limit oi g ^ 00. 
Equation (35) shows that the best extrapolation of the power is provided by the factor 
approximant, whose error is about 20%, while the root approximants have a twice larger 
error. The accuracy can be improved by defining the higher-order approximants. 

In the present section we demonstrated the application of the method of iterated root 
appproximants to the problem of calculating the ground-state energy of the anharmonic 
oscillator. This method works well for other problems too. For example, we have con- 
structed the iterated root approximants for the problem of the one-dimensional Bose gas 
of Sec. 3.2. The root approximant Rl{g), extrapolated from the left to right, gives the 
energy e(oo), as t — )■ 00, equal to 3.292, which is very close to the Tonks-Girardeau limit 
We have also considered the iterated root approximants for the problem of the 
diluted Fermi gas of Sec. 3.3. Extrapolating the ground-state energy from small g ^ 00 
to the strong-coupling limit — > 00 for Rl{g) provides an accuracy within the maximal 
error of order 10%. 

4.2 Odd factor approximants 

In the definition of the factor approximants of odd orders, there is a necessity of prescribing 
the value of one of the parameters, say Ai, in the general form (3). According to the scaling 
arguments [47,48] , this parameter can be set to one. Here we suggest one more variant of 
setting this parameter by defining it as ai/ao. 

Keeping in mind the same problem of extrapolating the ground-state energy of the 
anharmonic oscillator with Hamiltonian (26), we need to compare the accuracy of the 
corresponding factor approximants. The standard form is 

f^ig) = ^(1 + + ^2<7)"^(l + A,gr , (36) 
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where 

A2 = 21.86082, ^3 = 8.48018, 
ni = 0.27622, ria = 7.16531 x 10~^ = 0.12584 . 

In the strong-couphng hmit, this gives 

f*{g) ~ 0.669^°-4°9 ^ 00) . (37) 

Another variant of the factor approximant reads as 

frig) = ^ (1 + ^ (1 + A,gr-{1 + A,gr^ , (38) 
with the parameters 

A2 = 22.16875, ^3 = 8.83021, 
ni = 0.2212, 712 = 6.55016 x 10-^ rig = 0.11585. 

Now, the strong-couphng hmit becomes 

frig) = 0.718^°-^^^ (^ ^ 00) . (39) 

As is seen, the first form (36) extrapolates better the amphtude, while the second variant 
gives a better extrapolation of the power. Generally, there is no appriori preference for 
choosing this or that form, which, actually, is in agreement with the scaling arguments 
[47,48]. 



4.3 Weighted factor- root approximants 

One more possibility is to construct the weighted approximants defined as the linear 
combination of the factor and root approximants of the form 

Wtig) ^ Xflig) -r il - \)Rlig) . (40) 

with all parameters determined from the accuracy-through-order matching. The accuracy 
of such weighted approximants can be essentially improved. A more general way of 
constructing average values for an ensemble of approximants was considered in Ref. [66] 



5 Problem of extrapolation to infinity 

The self-similar factor approximants (3) extrapolate the asymptotic series (2), valid for 
small — 7> 0, to the region of finite (7 > 0. These approximants, as has been shown by 
numerous examples [45-49], provide for finite, even rather large, g very good approxima- 
tions, essentially more accurate than Pade approximants. 

But the problem, we face now, is to extrapolate series (2) not simply to finite or large 
g, but to find the limit /(oo) at — >■ 00, of course, assuming that this limit exists, so that 

/(oo) = lim f\g) = const. (41) 
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By requiring that the hmit 



/,*(oo)=lim/,*(^) = const (42) 



g^oo 



would also exist, in view of Eq. (3), we come to the condition 



limfoig)ll{AigP^ const. (43) 



1=1 



Let, for concreteness, the behavior of fo{g) at large g be fo{g) ~ Ag"' {g — >■ oo), as in 
Eq. (13). Then, for condition (43) to hold, it is necessary and sufficient that 



Nk 

a + Y^m^O. (44) 

i=l 

Therefore the value of approximant (3), given by 

fkioo) ^ AU AT , (45) 

i=l 

provides the approximation for the sought limit /(oo). 

This method of defining the limit f^{oo) at gf ^ oo by imposing the restriction (44) 
on the powers of the approximant (3) is very simple. However, as has been analyzed 
in Refs. [45,46], its accuracy is not high. In the following section, we suggest another 
method, whose high accuracy will be illustrated by calculating the pressure of fluctuating 
membranes. 



6 Method of variable transformation 

Instead of considering the limit /(oo) at g — )■ oo, it is convenient to make the change of 
variables 

9 = g{z) , z = z{g) , (46) 

such that 

Jim z{g) = 1 . (47) 

Then, for the function 

F{z) ^ f{g{z)) , (48) 

the sought limit is given by 

F(l) = /(oo) (49) 

at 2; = 1. 

With the change of variables (46), the series (2) become 

h{g{z))^Fk{z) (^^0), (50) 
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where 

k 

Fk(z) = Fo{z) bnz^ (51) 

ri-O 

is an expansion in powers of z up to the k-th order, with the coefficients bo = 1 and 
bn = bn{(ii, 02, ... , afc) defined through the coefficients a„ of series (2). Then, constructing 
for sum (51) the factor approximant F^{z), we obtain the approximation F^{1) for the 
sought hmit /(oo). 

In order to specify transformation (46), we assume the following natural properties. 
According to condition (47), it should be that z ^ 1 as g ^ oo. If the asymptotic 
behavior of the sought function is 

fig) ^ /(oo) (^1 + ^) (^ ^ oo) , (52) 

we require that transform (48) be 

F{z)c,F{l)[l + C2{l-z)] (z^l), (53) 

where Ci and C2 are constants and a; > 0. And, in agreement with Eqs. (50) and (51), 
we assume that z ^ as g ^ 0, so that 

g{z) c^Xz (z^O) , (54) 

with a scaling parameter A > 0. 

Comparing Eqs. (52) and (53) requires that 

§ c C2{1 -z) (^ ^ 1) , (55) 
9 

where we take into account that f{oo) = F{1). This yields 

where C3 = [Ci/ C2Y^^ ■ The interpolation between limits (54) and (56) can be done by 
using the self-similar factor approximants [45-50], which results in the form 

Thus, we obtain an explicit expression for the change of variables (46). For simplicity, we 
set in what follows the scaling parameter A = 1. 

But we need yet to define the exponent uj. For this purpose, we introduce the function 

_ d\Tif{g) _ g df{g) 

It is easy to notice that, if f{g) enjoys the asymptotic behavior (52), then the function 
(58) behaves as 

Pig) ~ - cuCig-'^ (g^oo). (59) 
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Therefore, the exponent u is defined as 

„ = (60) 

g->oo lag 

This procedure is to be accomphshed in each approximation order. That is, for the given 
fkig), defined in Eq. (2), we write 

dlnfkjg) , . 

^'^'^ ^ ^h^ ^''^ 

and expand this in powers of g, getting 

f3k{g) = f3o{g) E c„^" , (62) 

n=0 

with the coefficients Cq = 1 and c„ = c„(ai, 02, ... , a^) prescribed by the coefficients a„. 
Then, we construct the factor approximant 

Pl(g)^Po{g)U(^ + Digr' (63) 

i=l 

for series (62). Prom here, considering the hmit — > 00, we get 

.,._Un.!«, (64) 

in analogy with Eq. (60). 

Let, for example, the first factor in Eq. (63) behave as 

/3o{g) ~ Bg-^ (g^oo). (65) 

Then approximant (63), at large — > 00, is 

pk{g) ^ Bg^ n i^igT' (g^oo). (66) 

i=l 

As a result, Eq. (64) gives 

- ry + Ylrni\ ■ (67) 

To summarize, the calculational scheme is as follows. For a given series (2), we find 
the function (62) and construct its factor approximant (63), which provides us with the 
exponent (67): 

fkig) ^ Pkig) ^ Plig) ^ cu^ • (68) 

Then, making the change of the variable 
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from fk{g), we obtain expression (51). Constructing for the latter the factor approximant 
F^{z), we come to the value F^{1) approximating the sought limit /(oo): 

F,{z) ^ F^{z) ^ F;(1) . (70) 

When constructing the self-similar approximant F^{z), we need to solve the system 
of equations (9). Sometimes, though rare, it may happen that Eqs. (9), for some k, have 
no solutions which would yield real self-similar factor approximants. In that case, for 
Fj^{z) we take the arithmetic average of its neighbors F^ — {Fj^_^ + F^_,_^)/2. In practical 
calculations, one always deals with the approximation orders A; = 1, 2, . . . , X up to a finite 
maximal order K. Then the final answer for the set of F^ is given by the average of two 
last terms (FJ^ + F^_i)/2. The scheme, formulated in the present Section, will be applied 
to studying the properties of fluctuating membranes in the following Sections. 



7 Energy of one-dimensional membrane 

A cartoon of a membrane is a one-dimensional string oscillating between two rigid walls. 
This model, to our knowledge, was suggested by Edwards [67] and later considered in 
many articles [2-6,68]. It has been shown that calculating the free energy of the string is 
equivalent to finding the ground-state energy of a quantum particle in a one-dimensional 
box. Replacing the rigid walls by a soft potential, characterized by a finite stiffness g, 
and employing perturbation theory with respect to g yields the series 

2 k 

Ek{g) = ^,Y. ^ng-" (71) 

°y n=0 

for the particle ground-state energy E{g)^ with the coefficients 

1 1 1 1 

ao = 1 , ai = - , a2 = — , 03 



4 ' ^ 32 ' ^ 512 ' 

1 ^ 1 

04 = , 05 = -7—- , ag = , ay 



131072 ' o ' ' 16777216 
5 

(Zg — , flg — — 



8589934592 ' 
and so on. 

The scries (71) are obtained for the asymptotically small — )> 0. But, in order to pass 
to the case of hard walls, one has to consider the limit g ^ 00, with E{oo) being the 
sought value. Fortunately, the one-dimensional case allows for an explicit solution [67-69] 
giving 



from where 



^(00) = — = 0.077106 . (73) 
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This makes it possible to evaluate the accuracy of the self-similar approximants £^^(oo) 
with respect to the exact limit (73). 

We also wish to compare the accuracy of the two methods described above. First, we 
use the direct method by imposing restriction (44). To this end, we construct the factor 
approximants 

Eli9) = ^2^i^ + ^^9r (74) 
for series (71). Imposing the power restriction condition (44), we have 

E^i = 2. (75) 

i=l 

So that the sought limit is given by 

^^M = ^n^r, (76) 

^ i=i 

according to Eq. (45). The accuracy of approximants (76) is characterized by the per- 
centage errors 

^ Elioo) -E{oo) ^^^^^ 
^ E{oo) ^ ' 

with respect to the exact value (73). 

Another way is to follow the method of Sec. 6. Then we find the function (61), with 
expansion (62), which reads as 

k 

Pk{9)^-2f2cng\ (78) 



n=0 



with the coefficients 



_ 1 _ _ 1 



8' ' ' " 1024' 

3 5 

C4 = , C5 = - , C6 = , C7 = 



C8 = , Cg 

etc. 



262144 ' o ' ' 33554432 

35 



17179869184 ' 



Constructing the factor approximants (63) for series (78), we find 

2^ 



'1(9) = - 2 (79) 



V64 + g'' 

for all k > A. Therefore, the exponent (34) is 

Wjfc = 2 (A; > 4) . (80) 
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Hence, transformation (69) becomes 
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Using g = g{z), given by Eq. (81), we get 

F{z) ^ E{g{z)) , 
similarly to Eq. (48). And the related series (81) acquire the form 

^2 k 



TT 



n=0 



(81) 



(82) 



(83) 



with the coefficients 

bo = l 



64 = - 
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67 = - 



1024 ' 
77295 



3 

4 ' 

h = - 



62 = - 



32 ' 



= - 
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68 = - 



1185 
131072 ' 

119595 
33554432 ' 



&B = - 



16777216 ' 
and so on. 

Then, for series (83), we find the factor approximants 



512 ' 

1635 

262144 ' 
24489285 



8589934592 ' 



2 ATfe 



(84) 



whose values -^^(1) approximate the sought limit E{oo). The accuracy of F^{1) is char- 
acterized by the percentage errors 



e{F,) = ^(^) • 100% . 



(85) 



The results of our calculations for the factor approximants i5^(oo) and F^{1) are 
presented in Table 1, together with their errors (77) and (85). As is seen, the method of 
variable transformation of Sec. 6 is two orders more accurate than the method of power 
restriction employing restriction (44). The final answer, given by the former method, is 
E*{oo) = 0.0771, deviating only by 0.01% from the exact value (73). 



8 Pressure of fluctuating membrane 

An important class of membranes is formed by those membranes whose constituent 
molecules can freely move within them. Such membranes are called fluid. The ther- 
mal fluctuations of these membranes, at a temperature T, are controlled by their bending 
rigidity k. When modeling these membranes, one usually considers them as having a finite 
length L and an area A ^ 00. In order to describe their properties, one, first, assumes 
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that a membrane is located between the walls of a finite stiffness g. This allows one to 
resort to perturbation theory in powers of g. But to return to the case of hard walls, one 
needs to consider the limit g — ?■ oo, which requires to invoke a rcsummation procedure. 

It is convenient to introduce the dimensionlcss pressure p{g) of a fluctuating membrane, 
connected with the dimensional pressure P{g) through the relation 

Pi9) = ^ Pig) ■ (86) 
The asymptotic behavior of this function, at small gi ^ 0, is represented by the series 

2 k 

P^i9) = ^, Y.'^ng''. (87) 

°y n=0 

The coefficients of the perturbation series (87) are known only up to the sixth order [31], 
being 

ao = 1 , «i = ^ . " ^ ' " 2.176347 x 10"=^ , 



a4 



0.552721 X 10"^ , as = -0.721482 x 10"^ , ae = -1.777848 x 10"^ . 



We may notice that, up to the second order, the coefficients a„ in pressure (87) are the 
same as a„ in the ground-state energy (71). The pressure of the membrane, located 
between hard walls, is given by the limit 

p(oo) = lim p[g) . 

We shall again find this limit by two methods, by the method of the power restriction 
(44) and the method of the variable transformation of Sec. 6. 

In the direct method of power restriction (44), we find the factor approximants Pk{o6) 
corresponding to series (87). The approximants pg and pg cannot be defined by this way. 
And other approximants arc 

pt(oo) = 0.0193 , p;(oo) = 0.0232 , p*(oo) = 0.3120 , ^4(00) = 0.2880 . 

The most accurate Monte Carlo calculations for the membrane pressure have been ac- 
complished by Gompper and KroU [33] giving 

Pj^c = 0.0798 ± 0.0003 . (88) 

As we see, the accuracy of p^,(oo), compared to the Monte Carlo value (88), is rather bad. 

Now wc pass to the more elaborated method of Sec. 6, based on the change of variables 
prescribed by Eq. (69). The correct choice of the exponent uj^ is very important for 
achieving a good accuracy of the sought limit F^{1). This exponent is expressed by Eq. 
(64) through the function (63), for which we take the even factor approximants completely 
defined in Sec. 2. The series (62), for the considered case, has the same form as in Eq. 
(78), but with the coefficients 

Co = 1 , ci = - ^ , C2 = , C3 = 0.64173 x 10"^ , 
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C4 = 0.10668 X 10-^ , C5 = 0.46253 x 10"^ , cq = 0.18454 x 10"^ . 

Again, we may notice that the coefficients Cq, Ci, and C2 for the function (3k{g), in the case 
of the membrane, are the same as for series (78) in the case of the string. Constructing 
the even-order factor approximants f3l{g) and substituting these into Eq. (64), wc find 
that the sole real exponent cuk is given by the fourth-order approximant. Thus, we are 
left with the exponent 

Uk = 1.927 (A; > 4) . (89) 

Accomplishing the change of variables (69), we find series (51), for which we construct 
the factor approximants F^{z). Taking the limit 2; — > 1, we obtain 

F;(1) = 0.0906 , F5*(l) = 0.0898 , F6*(l) = 0.0747 . 

Averaging the last two values, we get our final result for the pressure (86) of the fluctuating 
membrane: 

p{oo) = 0.0823 . (90) 

This value is very close to the result [31] of Kastcning p{oo) = 0.0821, though it is 
3% higher than the Monte Carlo value (88) of Gompper and KroU [33]. The achieved 
accuracy is quite good, especially keeping in mind that the method of self-similar factor 
approximants is much simpler than the numerical method used by Kastening [31] and 
several orders simpler than the Monte Carlo simulations [33]. 

9 Conclusion 

We have suggested several modifications for constructing self-similar approximants in the 
frame of the self-similar approximation theory. Two main problems are considered, the 
problem of interpolation and extrapolation of asymptotic series. The suggested methods 
are illustrated by examples typical of chemical physics and quantum chemistry. 

A special attention is payed to the problem of defining the value of a function at infin- 
ity from its expansion at asymptotically small variables. We have designed a new way for 
constructing the self-similar factor approximants, so that to derive an accurate extrapo- 
lation /(oo) for a function f{g) in the limit of large g — > 00, when only the asymptotic 
series fk{g) at small g ^ are available. We have analyzed and compared two variants 
of the extrapolation. One of them involves a restriction on the powers of the constructed 
factor approximants, given by Eq.(44). This variant, however, is not sufficiently accurate. 
The latter is caused by the fact that the self-similar factor approximants are the most 
accurate when they are completely defined, by the re-expansion procedure, through the 
coefficients of the asymptotic series (2). But imposing additional constraints disturbs the 
self-consistency of the procedure and worsens the accuracy. 

The variant of Sec. 6, based on the variable transformation, is essentially more accu- 
rate. This is because it does not involve a restriction on powers. Vice versa, it takes into 
account the additional information on the behavior of the function f{g) when approach- 
ing the limit f{oo). The prescribed change of the variable, not merely tells that f{oo) is 
finite, but also describes how f{g) approaches this limit. The accuracy of the method is 
illustrated by calculating the pressure of fiuctuating fiuid strings and membranes. 
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In order to concisely summarize the ideas of the most accurate method, let us briefly 
delineate its main steps. Suppose we aim at flnding the limit /(oo) of a function f{g), 
as (? — 7> oo. But what is known for us is only the approximate behavior of the function 
at asymptotically small g ^ oo, where it is approximated by the series fk{g), and that 
f{g) — > const, as gr — > oo. For a given fkig), we define the function Pkig) through Eqs. 
(61) and (62). Then wc construct the factor approximants Pl{g), as in Eq. (63), and 
define the exponent cuk by Eq. (64). According to Eq. (69), we make the transformation 

^'^^^^ ^ (1 _ ^)iM ' 

and, as in Eq. (50), introduce Fk{z) = fk{gk{z)). Constructing the factor approximants 
F^{z) and taking the hmit z — )■ 1, we obtain the values F^{1) approximating the sought- 
function limit /(oo). Schematically, all this procedure is represented as the sequence of 
the following steps: 

fk{g) Pk{g) P*kig) ^ ^k ^ gk{z) 

^ Fk{z) ^ Fl{z) ^ Fl{l) ^ /(oo) . 

We have illustrated the above approach by calculating the pressure of fluid fluctuating 
membranes. The latter form a rather widespread important class of membranes studied in 
biology and chemistry [2-6,70]. The asymptotic series for the pressure were derived from 
Helfrich model [5]. The developed methods can be applied to other systems, where one 
needs to extrapolate the sought function from the region of asymptotically small variables 
to their flnite values. Moreover, the suggested methods make it even possible to flnd, with 
a good accuracy, the limit of the function at infinity. The advantage of the developed 
methods is their simplicity and high accuracy. 
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Table 1 



k 


Elioo) 


% 




% 


3 


0.1500 


94.1 


0.0593 


-23.1 


4 


0.1370 


78.1 


0.0935 


21.3 


5 


0.0526 


-31.8 


0.0926 


20.0 


6 


0.0550 


-28.7 


0.0829 


7.52 


7 


0.1030 


33.4 


0.0732 


-5.01 


8 


0.0993 


28.8 


0.0805 


4.36 


9 


0.0620 


-19.6 


0.0803 


4.17 


10 


0.0636 


-17.5 


0.0783 


1.49 


11 


0.0926 


20.1 


0.0762 


-1.24 


12 


0.0906 


17.5 


0.0780 


1.13 


13 


0.0662 


-14.1 


0.0779 


1.08 


14 


0.0674 


-12.6 


0.0774 


0.37 


15 


0.0882 


14.4 


0.0768 


-0.34 


E*{oo) 


0.0778 


0.91 


0.0771 


-0.008 



The factor approximants El{oo) and -^^(1), together with their percentage errors, ap- 
proximating the ground-state energy E{oo). 
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